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Abstract 

We consider long-range self-avoiding walk, percolation and the Ising model on % d 
that are defined by power-law decaying pair potentials of the form D(x) x |x|~ rf ~ a 
with a > 0. The upper-critical dimension d c is 2(aA2) for self-avoiding walk and the 
Ising model, and 3(a A 2) for percolation. Let a ^ 2 and assume certain heat-kernel 
bounds on the n-step distribution of the underlying random walk. We prove that, 
for d > d c (and the spread-out parameter sufficiently large), the critical two-point 
function G Pc (x) for each model is asymptotically C\x\ aA2 ~ d , where the constant 
C € (0, oo) is expressed in terms of the model-dependent lace-expansion coefficients 
and exhibits crossover between a < 2 and a > 2. We also provide a class of random 
walks that satisfy those heat-kernel bounds. 
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1 Introduction 



The two-point function is one of the key observables to understand phase transitions and 
critical behavior. For example, the two-point function for the Ising model indicates how 
likely the spins located at those two sites point in the same direction. If it decays fast 
enough to be summable, then there is no macroscopic order. The summability of the 
two-point function is lost as soon as the model parameter (e.g., temperature) is above the 
critical point, and therefore it is naturally hard to investigate critical behavior. 

The lace expansion is a powerful tool to rigorously prove mean-field behavior above 
the model-dependent critical dimension. The mean-field behavior here is for the two-point 
function at the critical point to exhibit similar behavior to the underlying random walk. 
It has been successful to prove such behavior for various statistical-mechanical models, 
such as self- avoiding walk, percolation, lattice trees/animals and the Ising model. The 
best lace-expansion result obtained so far is to identify an asymptotic expression (= the 
Newtonian potential times a model-dependent constant) of the critical two-point function 
for finite-range models, such as the nearest-neighbor model. However, this ultimate goal 
has not been achieved, before this paper, for long-range models, especially when the 1- 
step distribution for the underlying random walk decays in powers of distance; only the 
infrared bound on the Fourier transform of the two-point function was available. This 
was partly because of our poor understanding of the long-range models in the x-space, 
not in the Fourier space. For example, the random-walk Green's function is known to be 
asymptotically Newtonian/Riesz depending on the power of the aforementioned power- 
law decaying 1-step distribution, but we were unable to find optimal error estimates in 
the literature. Also, the subcritical two-point function is known to decay exponentially 
for the finite-range models, but this is not the case for the power-law decaying long-range 
models; as is shown in this paper, the decay rate of the subcritical two-point function is 
the same as the 1-step distribution of the underlying random walk. 

Therefore, the goal of this paper is to overcome those difficulties and derive an asymp- 
totic expression of the critical two-point function for the power-law decaying long-range 
models above the critical dimension, using the lace expansion. We would also like to in- 
vestigate crossover in the asymptotic expression when the power of the 1-step distribution 
of the underlying random walk changes. 



1.1 Models and known results 

Self-avoiding walk (SAW) is a model for linear polymers. We define the two-point function 
for SAW on Z rf as 



\u 

g s p aw (x) = p h n d ^ - w i-o ri( i - ^ 

lo:o— >X j=l 8<t 



where p > is the fugacity, \u\ is the length of a path uj = (ujq,uji, . . . ,00^) and D : 
Z d — > [0, 1] is the Z d -symmetric non-degenerate (i.e., D(o) ^ 1) 1-step distribution for 
the underlying random walk (RW); the contribution from the 0-step walk is considered 
to be 5 QjX by convention. If the indicator function T^f s<t (1 — S Ul>tUt ) is replaced by 1, 
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then Gp AW (x) turns into the RW Green's function G p w (x), whose radius of convergence 

is 1, as X p w = ExezJ G p W ( x ) = i 1 ~ P)' 1 for P < 1 and X™ = oo for p > 1. 
Therefore, the radius of convergence p^ AW for G p AW (x) is not less than 1. It is known that 
Xp AW = ^2 xe z d Gp AW (x) < oo if and only if p < p^ AW and diverges as p j" p^ Aw . Here, and 
in the remainder of the paper, we often use "=" for definition. 

Percolation is a model for random media. Each bond {u, v}, which is a pair of vertices 
in Z d , is either occupied or vacant independently of the other bonds. The probability that 
{u, v} is occupied is defined to be pD{y — u), where p > is the percolation parameter. 
Since D is a probability distribution, the expected number of occupied bonds per vertex 
equals pE^o^W = p(1 ~~ ^(°))- The percolation two-point function G p mic (x) is defined 
to be the probability that there is a self-avoiding path of occupied bonds from o to x; 
again by convention, G£ erc (o) = 1. 

The Ising model is a model for magnets. For Ac Z d and <p = {y9„}„eA £ {il}^ we 
define the Hamiltonian (under the free-boundary condition) as 

H\{v) = - j u ,v<fu<fv, (i-2) 

(u,d}cA 

where J U|l) = J a ,v-u > is the ferromagnetic pair potential and inherits the properties 
of the given D, as explained below. The finite- volume two-point function at the inverse 
temperature (3 > is defined as 

v>e{±i} A / ^6{±i} A 

It is known that (<p i Px)i3,A is increasing in A t Z d . Let p = XLez d tanh(/3J 0jX ). The Ising 
two-point function Gp ing (a;) is defined to be the increasing-volume limit of ((p <Px)p,A' 

G l ; ins (x) = hm(ip (p x ) PA . (1.4) 

' A^Z d 

Let -D(x) = p _1 tanh(/3 J 0jX ). 

For percolation and the Ising model, there is a model-dependent critical point p c > 1 
(from now on, we omit the superscript, unless it causes any confusion) such that 

t&A l = °° IP^Pc]' V |a,Ko0 l >0 \P>Pc\- 

The order parameter 6^ crc is the probability that the occupied cluster of the origin is 
unbounded, while 8 p ins is the spontaneous magnetization, which is the infinite- volume limit 
of the finite- volume single-spin expectation ((p )p A under the plus-boundary condition. 
The continuity of 9 P at p = p c in a general setting is still a remaining issue. 

We are interested in asymptotic behavior of G Pc (x) as \x\ — > oo. For the "uniformly 
spread-out" finite-range models, e.g., D{x) = l{| z | = i}/(2d) or D(x) = 1^\ x \\ oc <l}/(2L + l) d 
for some L G [1, oo), it has been proved [TBI E2 EE] that, if d > 4 for SAW and the Ising 
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model and d > 6 for percolation, and if d or L is sufficiently large (depending on the 
models), then there is a model-dependent constant A (= 1 for RW) such that 



(L6) 



where "~" means that the asymptotic ratio of the left-hand side to the right-hand side is 
1, and 



dT{^ 



27r d/2 ' 



a 2 = \A 2 D(x) = 0(L 2 ). (1.7) 



This is a sufficient condition for the following mean-field behavior (TJ El [31 El EI] : 

^ f \-i a ~ fVP-Pc [king], , s 

X P ~ (Pc -p) , #p ~ < r . ,. , (1-8) 

ptpc pJ-Pc I p — p c [percolation] , 

where "x" means that the asymptotic ratio of the left-hand side to the right-hand side is 
bounded away from zero and infinity. 

The proof of the above result is based on the lace expansion (e.g., [IHl EU |26] ) . The core 
concept of the lace expansion is to systematically isolate interaction among individuals 
(e.g., mutual avoidance between distinct vertices for SAW or between distinct occupied 
pivotal bonds for percolation) and derive macroscopic recursive structure that yields the 
random- walk like behavior ( II. 6p . When d > d c and rfVL > 1 (i.e., d or L sufficiently 
large depending on the models), there is enough room for those individuals to be away 
from each other, and the lace expansion converges [191 EH EE] . The resultant recursion 
equation for G p is the following: 

^pD(v)G p (x-v) [RW], 

v&Z d 

7Tp(x) + ix p (u) pD(v — u) G p {x — v) [Ising & percolation], 

u,vez d 

where tt p is the lace-expansion coefficient. To treat all models simultaneously, we introduce 
the notation / * g to denote the convolution of functions / and g in Z d : 

(f*g){x) = Y,f(v)g(x-v). (1.10) 

Then the above identities can be simplified as (the spatial variables are omitted) 

(5 + pD*G p [RW], 
G p = I 5 + (pD + n p ) * G p [SAW], (1.11) 

^7r p + Tip * p(D — D(o)8) * G p [Ising & percolation]. 



G p (x) 



4 



Repeated use of these identities yieldsH 

g p = n p + n p * pD * g p , 

where 

(S , x [RW], 

oo oo 



;i.i2) 



n p {x) 



n=0 

oo 



n=0 



n-fold 



;i.i3) 



^(-pD(o)r\; n (x) 



[Ising & percolation], 



^ n=l 



with the convention f*°(x) = 5 ^ x for general /. When d > d c and d\J L ^> 1, there 
is a p > such that \II Pc (x)\ is summable and decays as \x\~ d ~ 2 ~ p [161 HI I2S]. The 
multiplicative constant A in (ll.6p and p c can be represented in terms of II Pc (x) as 



Pc 



Pc 



x, 



;i.i4) 



In this paper, we investigate long-range SAW, percolation and the Ising model on Z d 
defined by power-law decaying pair potentials of the form D(x) x with a > 0. 

For example, as in [10], we can consider the following uniformly spread-out long-range D 
with parameter L G [1, oo): 



D(x) 



x_ 111 —d—a 

zJli 



E 



y_ in — d— a ' 
illll 



where II x|L = \x\ V I. As a result, 



£>(x) = 0(L a )|||x|| 



— d— a 
L i 



;i.i5) 



;i.i6) 



1 For SAW, since ||7T p ||i = o(l) asdVl^oo and ||G p ||oo < oo for every p < p c [16j ITS]. 
G p = 5 + pD * G p + n p * G p = 5 + pD * G p + n p * (5 + pD * G p + n p * G p ) 



replace 



{5 + 7T P ) + (5 + 7T P ) *pD*G p + n* 2 * Gp = 

replace 



(fTT2l) . 



For percolation and the Ising model, since D(o) = o(l) and p||7r p ||i = 1 + o(l) as d V L — > oo and 
HGpHoo < 1 for every p<p c [THHKIB], 

G p = 7t p + 7T P * p_D * Gp — pD(o)ir p * Gp 

replace 

= 7Tp + 7Tp * * G p — pD{o)-K p * {tt p + TTp * pD * G p — pD(o)n p * Gp) 

= (tt p - P D(o)t:; 2 ) + (tt p - pD(o)tt* p 2 ) *pD*G p + (- P D(o)) 2 tt* p 2 *G^ = ► (ffHa)- 

replace 
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which we require throughout the paper (cf., Assumption [T7T] below). The goal is to see 
how the asymptotic expression (jl.6p of G Pc (x) changes depending on the value of a. We 
note that (II. 6p and (I1.14p are invalid for a < 2 because then a 2 = oo. 
Let 

d = (2(aA2) [SAW rising], 
1 3 (a A 2) [percolation]. 

It has been proved [20] that, for d > d c and L ^> 1, the Fourier transform G p (k) = 
Ylx<=z d e lk x G p (x) for the long-range models is bounded above and below by a multiple of 
G^i^k) = (1 — pD{k))~ l with p = p/p c , uniformly in p < p c . Although this gives an 
impression of the similarity between G Pc (x) and Gf w (x), it is still too weak to identify the 
asymptotic expression of G Pc (x). The proof of the above Fourier-space result makes use 
of the following properties of D that we make use of here as well: there are v a = 0(L aA2 ) 
and e > such that 

Si ll°grH + 0(i) [" = 2], 

If a > 2, then f Q = a 2 / (2c/). Moreover, if L 1, then there is a constant A G (0, 1) such 
that| 

n<0(L-«)n-^ [„>1], 1 -/>(*)(< 2 " A ^ [ 7V^' ( h2 °» 

All those properties hold for D in (05]) (cf., [rO|fT2]). 
1.2 Main result 

In addition to the above properties, the n-step transition probability obeys the following 
bound: 



D *n (x)< 9i^l nx 1 [^ 2 ]' (1 2i) 

( } INir 2 iloglNU [a = 2]. 



This is due to the following two facts: (i) the contribution from the walks that have at 
least one step which is longer than c|||x|||l for a given c > is bounded by 0(L q )?t,/||x||^ +q ; 
(ii) the contribution from the walks whose n steps are all shorter than c|||x|||l is bounded, 
due to the local CLT, by 0(vn)~ d / 2 e-^ 2/ ° {i}n) < 0(vn)/l\x\\ d L +2 (times an exponentially 



2 In the proof of the bound on ||-D* n ||oo, we simply bounded the factor log ^ in PHI (A. 4)] by a positive 
constant. If we make the most of that factor instead, we can readily improve the bound for a = 2 as 

P* n ||oo <0(L- d )(nlogn)- d / 2 . (1.19) 
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small normalization constant), where v is the variance of the truncated 1-step distribution 
D(y) = ^d/)l{|j,|<c|x|} an d equals 

(\\x\\ 2 L - a [a<2], 
v = \v\ 2 D{y) < 0(L aA2 ) x I log \\x\\ L [a = 2], (1.22) 
y^ d [l [a>2\. 

For a^2, the inequality ( 11.21 jl is a discrete space-time version of the heat-kernel bound 
on the transition density p s (x) of an a-stable/Gaussian process: 

In Section |2~T| we will show that the properties ( ll.lfip . (jl.l8p and (ll.2ip are sufficient 
to obtain an asymptotic expression of Gf w (x). However, these properties are not good 
enough to fully control error terms arising from convolutions of D* n (x) and n p (x) in (11.131) . 
To overcome this difficulty, we assume the following bound on the discrete derivative of 
the n-step transition probability: 



D* n (x) 



D* n (x + y) + D* n (x - y) 



< ^2SK [u<m- d-24) 



Fill I, 



Here is the summary of the properties of D that we use throughout the paper. 

Assumption 1.1. The 7* d -symmetric 1-step distribution D satisfies the properties U.16]) . 
/ f05|) . CT . (EH]) and [T24\ ). 

In Appendix |AJ we will show that the following D satisfies all properties in the above 
assumption: 

D(x) = J2U! t (x)T a (t), (1.25) 

teN 

where Ul is in a class of Z d -symmetric distributions on Z d fl [-L, L] d , and T a is the stable 
distribution on N with parameter a/2 ^ 1. 

Under the above assumption on we can prove the following theorem: 

Theorem 1.2. Let a > 0, a ^ 2 and 

r(^p) 

2aA2 7r <i/2p^aA2 



7. = ^_J^L 2V (1-26) 



and assume all properties of D in Assumption Then, for RW with d > a A 2 and 
any L > 1, and for SAW, percolation and the Ising model with d > d c and L>1, there 
are fi £ (0,a A 2) and A = A(a,d,L) £ (0, oo) (A = 1 for random walk) such that, as 
\x\ — > oo, 

P A\x\ d - aA2 \ x \d-aA2+n ■ K 1 ^') 



7 



As a result, by IWtf . \ p and 9 P exhibit the mean-field behavior U.8\) . Moreover, p c and A 
can be expressed in term of II P in U.13\) as 



[a < 2], 

P ^n P M-\ A = * + <g2>|'iU.) [a>2]. (L28) 

X 

Remark 1.3. (a) The finite-range models are formally considered as the a = oo model. 
Indeed, the leading term in (I1.27P for a > 2 is identical to (II. 6p . 

(b) Following the argument in [TBI |2B] , we can "almost" prove Theorem 11.21 for a > 2 
without assuming the bounds on D* n (x). The shortcoming is the restriction d > 10, 
not d > 6, for percolation. This is due to the peculiar diagrammatic estimate in [16], 
which we do not use in this paper. 

(c) The asymptotic behavior of G Pc (x) in ( II. 6p or ( I1.27P is a key element for the so-called 
1-arm exponent to take on its mean-field value [T7 ] |2T | [23 ] 125] . For finite-range critical 
percolation, for example, the probability that o G Z d is connected to the surface of 
the d-dimensional ball of radius r centered at o is bounded above and below by a 
multiple of r -2 in high dimensions [22] • The value of the exponent may change in a 
peculiar way depending on the value of a 



(d) As described in (ll.28p . the constant A exhibits crossover between a < 2 and a > 2; in 
particular, A = p c for a < 2 (cf., f)3.6p below). According to some rough computations, 
it seems that the asymptotic expression of G Pc (x) for a = 2 is a mixture of those for 
a < 2 and a > 2, with a logarithmic correction: 



G ^ , , 1 | r (1-29) 

|x|-k» p c \x\ a z log \X\ 



One of the obstacles to prove this conjecture is a lack of good control on convolutions 
of the RW Green's function and the lace-expansion coefficients for a = 2. As hinted 
in the above expression, we may have to deal with logarithmic factors more actively 
than ever. We are currently working in this direction. 

1.3 Notation and the organization 

From now on, we distinguish G p w from G p for the other three models, and define 

S p = G p w . (1.30) 

Here, and in the remainder of the paper, the spatial variables are sometimes omitted. For 
example, 

S p = 5 + pD * S p (1.31) 
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is the abbreviated version of the convolution equation 

s P (x) = 5 , x + ( p d * s p )(x) = s 0>x + £ P D (y) s p( x - v)- (i- 32 ) 

y&Z d 

We also recall the notation 

= |z| VI. (1.33) 

The remainder of the paper is organized as follows. In Section El we prove the asymp- 
totic expression (11.271) for Si, as well as bounds on S p for p < 1 and some basic properties 
of G p for p < p c . Then, by using these facts and the diagrammatic bounds on the lace- 
expansion coefficients in [181 EE] , we prove d 1 . 2 7[) for G Pc in Section [3j 



2 Preliminaries 

In this section, we derive the asymptotic expression (I1.27P for Si, which will be restated as 
Proposition 12.11 an d prove some properties of G p that will be used to prove Theorem 11.21 
in Sectional 



2.1 Asymptotics of S p 

Proposition 2.1. Let q > 0, a ^ 2 and d > a A 2, and assume all properties but &1-24 ) 
in Assumption Then there is a fi G (0, a A 2) such that, for any L > 1, p < 1 and 

K > 0, 

6 , x <S p (x)<S 0>x + \ d _ a J [VxeZ d ], (2.1) 

I x I\l 



*(*) = ^ + ^ZZZ ] N>i' + 1. (2-2) 
where a constant in the 0(L~ aA2+M ) term depends on n. 

Proof. The inequality (12.11) is an immediate result of (OTP , p < 1 and fPUD - (TPTD aa^l 

0<S p (x)-5 o , x <J2D* n ^)<^^ E n + 0{L- d ) £ n~"W 

n=l \\ X \\l n= i n=(\\x\\ L /L)<** 2 

Q(^° A2 ) , 23) 

\\\X T 



3 For a = 2, we can readily bound S p (x) — 5 . x by using p,19[) for n > N x = \\x\W/(L 2 log |||x|||l) and 
([TT2TT) for n < N x as 

SJx) -5 ox < V D* n (x) + V D* n (x) < '- . 

~1 „fe. \lx\H- 2 \0gl\xl\ L 
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To prove the asymptotic expression (12. 2p . we first rewrite Si(x) for d > a A 2 as 
Sl ( x ) = [ e ~T = Tdt [ -^L e -**-*(i-A(*)) 



/ d t / A*"*"*" t(l-D(k)) dt + Iu (24) 



for any T G (0, oo), where 

y {2ny y ^ n! ^ ^ 

Next we rewrite the large-t integral as 

poo p jdi _ poo 5 

/ dt/ AJL-^-»m = dW*) + E J ;> (2-6) 



i=2 



where p s (x) is the transition density of an a-stable/Gaussian process (cf., (jl.23p ). and for 
any R G (0, 7r), 

h = - f T dt Pvat (x) = ~ [ T dt f -^-e-^-^\ (2.7) 

JO JO JR d \ £1X ) 

J 3 = Tdt [ -f^ e~ ik - x ( e -*( 1 -^( fe )) - e -M*r Aa> ) , (2.8) 



(2vr) d 



T ./|fc|<R 



T 



^4=/ dt/ 7 _ 3;fe -**-*a-oW)i {|ifc|>J9 , (2.9) 



[-7r,7r]d (2vr) c 

qA2 

( O \d 

T J\k\>R 



d d k 



h = - I dt I ^_ e -ik-x-v a t\k\^\ ( 210 ) 



dte-^l^ = = - J^ / dt t^V^, (2.11) 



By using the identity 

-v a t\k\ aA2 _ 

' Va \k\^ v a r(^) y 

we obtain 

oo 

qA2 







u«r(^) y " y^ (27r) d _ |x| 



ttvM = -^tW dt t-- f - e-l*^ = 1 J&. (2.12) 



As a result, we arrive at 

5 

a/ i 
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It remains to estimate Ej=i J r First > h Y CM) and dH2S]) , we can estimate Ii + J 2 
for I a? I > L as 



0(L aA2 ) f T J 0{L aA2 )T 2 



l / i + / 2 |<Tn^/ dtt< ; 2 . (2.14) 



|.r| 
Let 

,.JMl T _fMV A "" ! (215) 

Then, we obtain 

| /l + /2 l<£(^). (2 . 16) 

Next we estimate 1%. For small R, whose value will be determined shortly, we use 
([USD to obtain 



e -t(l-U(fc)) _ e -v a t\k\^ 2 



< 0(L aA2+e )t\k\ aA2+e e- v « maA2 . (2.17) 



Therefore, by (12.151) . 

poo /' 

|/ 3 1 < 0(L QA2+e ) / dt t / d d k \k\ aA2+ * e - Vam 

JT J\k\<R 



\k\<R 

poo 

0(L aA2+e ) / dt t I — ( — ] e 



on pd tR aA2 1 / \ d + aA2 + e 

dr 1 r \ aA2 



T JO 



r \v a t 



<o{l« a2+ <) I dtt{ Va ty i± ^ 

T 

< 0(L-)T-Si = (2.18) 



Finally we estimate J4 + I5 and determine the value of R during the course. First, by 
(P8|) - ffl~20]) . we have 



d d k 



\h\< J t dtj ^^^e-^^lOT^^^pii^^^j + ^iifc^^-!}) 

£ r d KL(&^ +o(iK,a ) 

POO 

<0(L- d ) / dtt-^T(^;tc(LR) aA2 ) + 0(l)e- TA , 
Jt 



(2.19) 



where T(a;x) = dtt a 1 e 1 is the incomplete gamma function, which is bounded by 
0{x a ~ 1 )e~ x for large x. Here we choose R to satisfy 

tc(LR) aA2 = log t. (2.20) 

a A 2 



11 



Then, for large t, 

r(aA2' tc(LR) aA2 ) < 0((tc(LR) aA2 )^- l )e- tc ^ aA2 

= 0{(\ogt)^- l )r^ < 0(t"&3). (2.21) 

Therefore, again by ( )2T5]) (cf., ( |2T8|l ). 

, dt HLRr 2 ) < 0{L- d )T^ = ^_ aA2 J . (2.22) 

We can estimate J5 in exactly the same way. The exponentially decaying term in 
(I2.19P obeys the same bound, since, for sufficiently large N (depending on k), 

„—TA ^ CN y-drpl-^ T drp~(N+l— ^±|) 

< c N L~ d T l ~^L d ~ {N+l ^^ ){ah2 ~^ /2)K (■; \x\ > L 1+K T > L (aA2 "^ /2)K ) 



< c N L~ d T^ = ^— —L . (2.23) 



Summarizing the above, we obtain that, for |x| > L 1+K , 
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< — -. (2.24) 



3=1 

This together with ( 12 . 1 3[) completes the proof of Proposition 12.11 



2.2 Basic properties of G p 

In this subsection, we summarize some basic properties of G p . Roughly speaking, those 
properties are the continuity up to p = p c (Lemma 12. 2p . the RW bound that is optimal 
for p < 1 (Lemma I2.3P and the a priori bound that is not sharp but finite as long as 
p < p c (Lemma 12. 4p . We will use them in the next section (especially in Section [372]) to 
prove Theorem 11.21 

Lemma 2.2. For every x G 'L d , G p (x) is nondecreasing and continuous in p < p c for 
SAW, and in p < p c for percolation and the Ising model. The continuity up to p = p^ AW 
for SAW is also valid if G SAW (x) is uniformly bounded in p < p^ AW . 

Proof. For SAW, since G p AW (x) is a power series of p > with nonnegative coefficients, 
it is nondecreasing and continuous in p < pl Aw . The continuity up to p = p^ AW under the 
hypothesis is due to monotone convergence. 

For the Ising model, we first note that, by the Griffiths inequality [13], {(p <£x)i3,A 
is nondecreasing and continuous in > and nondecreasing in A C Z d . Therefore, 
the infinite-volume limit G p irLS (x) is nondecreasing and left-continuous in p > 0. The 
continuity in p < p l * ins follows from the fact that, for p < p^ [ng , G l p sins (x) coincides with 
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the decreasing limit of the finite-volume two-point function under the "plus-boundary" 
condition, which is right-continuous in p > 0. 

For percolation, G^{x) is nondecreasing in p > because the event that there is a 
path of occupied bonds from o to x is an increasing event. The continuity in p > is 
obtained by following the same strategy as explained above for the Ising model and using 
the fact that there is at most one infinite occupied cluster for all p > 0. This completes 
the proof of Lemma 12.21 ■ 

Lemma 2.3. For every p < p c and x G Z d , 

G p {x) < S p {x), P D(x)(l - 8 , x ) < G p (x) - 5 Q)X < (pD * G p )(x). (2.25) 

Proof. The first inequality for p > 1 = p™ is trivial since S p (x) = oo for every x G Z d . 
On the other hand, the first inequality for p < 1 is obtained by using the second inequality 
N times and then using (11.201) . as 

N-l 

G p (x)<Y^p n D* n (x)+p N (D* N *G p )(x)<S p (x) + \\D* N \\ ooXp -> S p (x). (2.26) 

n=0 

It remains to prove the second inequality in ( 12.251) . In fact, it suffices to prove the 
inequality only for i^o, since G p [o) = 1 for all three models and therefore the inequality 
is trivial for x = o. For SAW and percolation, the inequality is obtained by specifying the 
first step pD and then using subadditivity for SAW or the BK inequality for percolation 
[6] . For the Ising model, we use the following random-current representation [lj [TJ] (see 
also [261 Section 2.1]): 

dn={o}A{x} -i-r {@Ju,v) nu,V fr> „n 

WoVx)p,A = ^ — — , w A {n) = ■ — , (2.27) 

dn=0 

where n = {n u>v } is a collection of Z + -valued undirected bond variables (i.e., n U)V = n VjU G 
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>).-- 



Z + = {0} U N for each bond {u,v} C A), dn is the set of vertices y such that 5^ 2eA 
is an odd number, and "A" represents symmetric difference (i.e., {o}A{x} = if x — o, 
otherwise {o}A{x} = {o,x}). Using this representation, we prove below that, for i^o, 

pD(x) < (<Po¥x)p,K < S ^pD(y){^ x )p A , (2.28) 

yeA 

where pD(x) = tanh(/3J D x ). The second inequality in (I2.25P for the Ising model is the 
infinite- volume limit of the above inequality. 

To prove the lower bound of (I2.28p . we first specify the parity of n ox to obtain that, 
for x 7^ o (so that {o}A{x} = {o,x}), 

Y w A (n) + Y w a(^) 

dn={o,x} dn={o,x} 
, v ( n o,x odd) (n 0tX even) 

{<Po(Px)PA = 7— ^ — ■ ( 2 - 29 ) 

> J WA{n) + w A{n) 

dn=0 dn=0 
(n ,x odd) (n ,x even) 
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Let 

Y y (z,x)= £ w A(n), Z y = £ w A(n). (2.30) 

dn={z}A{x} dn=0 
(n , y even) ( n o,y even) 

Then, by changing the parity of n ox (and the constraint on On accordingly) and recalling 
tanh(/3J oa .) — pD(x), we obtain 

£ w A (n)= P D(x) J2 VA(n)=pD(x)Z x , (2.31) 

dn={o,x} dn=0 
(n , x odd) (n ,x even) 

w A (n)= P D{x) Yl w A (n) =pD{x)Y x {o,x), (2.32) 



dn=0 dn={o,x} 
(na, x odd) (n o x even ) 



hence 



pD(x)Z x +Y x (o,x) , (1 - p 2 D(x) 2 )Y x (o, x) . 

(<Po<Px)f),A = ~ = p£>(x) + ^ ~ L > pD(z). (2.33) 

p£)(x)y r (o, x) + Z x pD(x)Y x (o, x) + Z x 

To prove the upper bound in f)2.28p . we first note that, if dn = {o, x}, then there must 
be at least one y G A such that n 0)y is an odd number. By similar computation to f)2.3ip . 
we obtain that, for i^o, 

Mn)<J2 E Mn) = J2pD(y) £ w A (n) . (2.34) 

dn={o,x} y&A dn={o,x} y&A dn={y}A{x} 

(n 0i y odd) ("o,;/ even) 

11 v ' 



On the other hand, J2d n =0 w A(n) > ^» f° r an y 2/ £ A. Therefore, for i^o, 
£ ™A(n) 

- ^ ; <£ ^p'^ < £^(,)(^ )/3 , A , (2.35) 
2^ w a(«) <,eA Z ^ ye A 

dn=0 

where we have used the fact that Y y (y,x)/ Z y is equivalent to the finite-volume two-point 
function under the restriction J 0>2/ = and therefore, by using the Griffiths inequality [13], 
it is bounded above by ((p y (p x ) 0,a- This completes the proof of (12.281) . hence the proof of 
Lemma 12.31 ■ 

Lemma 2.4. Assume the property U.16]) in Assumption Then, for every a > and 
p < p c , there is a K p = K p (a, d, L) < oo such that, for any x G 7L d , 

G p (x) < K p \lxl\- L d - a . (2.36) 

Remark 2.5. This together with the lower bound in (I2.25P implies that, for every p < p c , 
G p (x) is bounded above and below by a p-dependent multiple of |||x|||^ d_a . This shows 
sharp contrast to the exponential decay of G p {x) for the finite-range models. 
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Proof of Lemma \KQ Since G p (o) < x P < oo for p < p c , it suffices to prove (I2.36P for 
i^o. We follow the idea of the proof of [H Lemma 5.2] for one- dimensional long-range 
percolation and extend it to those three models in general dimensions. The key ingredient 
is the following Simon-Lieb type inequality: for < t < \x\, 

G p (x)< G p (u)pD(v-u)G p (x-v). (2.37) 

{u,v}cZ d 

(M<«M) 

For SAW and percolation, this is a result of subadditivity or the BK inequality (cf., e.g., 
[T5"| 124"] ) . For the Ising model, this is obtained by using the random- current representation 
( I2.27P and a restricted version of the source-switching lemma (26J Lemma 2.3], as follows. 
Let Z\ = J2dn=0 w ^( n ) sucn that, for i^o, 

<y^W= £ (2.38) 

dn={o,x} 

We note that, if dn = {o, x}, then there is a path u = (uq, u\, . . . , Ut) C A from ojq = o 
to oj t = x such that n U)s _ ltbJs is odd for every s £ {1, . . . ,t}; moreover, there is a unique 
r G {1, . . . ,t} such that |u> r _i| < I < \u T \ (i.e., r is the first time when u crosses the 
surface of the ball Bi of radius I centered at the origin). This can be restated as follows: if 
dn = {o, x}, then there is a bond {u, v} C A such that n Uj „ is odd and that u is connected 
from o with a path of bonds C B^ with odd numbers. Therefore, 

/ \ -v x -v w\(n) 

VPo<Px)P,K< y. }^ Mn u , v odd} l{o^->u in B e }, (2.39) 

z — ' z — ' Z/\ « 

{«,»}CA 9n={o,i} 

(\u\<e<\v\) 

where {o <; — > u in S^} is the event that o is connected to u with a path of bonds b C Bg 

n 

satisfying n& > 0. Multiplying Zr^/Zr^ = 1 to both sides of (I2.39P and using the identity 
z s t = Y.dm=z w B l { m ), we obtain 

hn,n\ <r \" w B.(m) w A (n) 

VPoVx)P,h S > > ^ ^ l{n„,„ odd} ^{o^u in S^}' l^- 4U J 

Z/R. Z/A m+n 

{ii,u}CA 3m=0 ^ A 

(\u\<£<\v\) dn={o,x} 

where we have used the trivial inequality t{ oi m iaB t } < l{o< >Min_B f }- Then, by using 

n m-\-n 

the source-switching lemma [261 Lemma 2.3], we obtain 
/ \ / \^ w Bl (m) w A (ri) 

\</V^//3,A S ^ 7 JL{n„,„odd} JL{o^uinBa 

Z R. Z/A rra + rx 

{u,ii}cA am={o}A{u} £ 
(M<^<M) <9ra={u,2;} 

= £ w«>a* £ ^> ( 2 - 41 ) 

{■u,ij}CA cto.={u,:c} 
(M<*<M) (%,, odd) 
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where we have used the identity 1{ < >u in B( y = 1 given dm = {o}A{u} and then used 

TTl + Tl 

(12.270 . Finally, by following the same argument as in (I2.34p -f l2.35p and then taking the 
infinite- volume limit, we obtain (I2.37P for the Ising model. 

Now we prove (12.360 by using (12.371) with £ = (the factor | is unimportant as 
long as it is less than |). Let 

c x = G p (u)D(v-u). (2.42) 



{u,v} 

(M<i|x|<M) 

We note that c x — > as |a;| — > oo, because 

c x = G p (u)pD(v — u) + G p (u)pD(v — u) 

{u,v}dZ d {u,v}dZ d 
(\u\<l\x\, §M<M) (±\x\<\u\<i\x\<\v\) 

<X P P sup D(v-u)+p G p( u )- ( 2 - 43 ) 

U '\ U \-i\ X \ v.\v\>±\x\ u:\u\>\\x\ 

" v ' v v ' 

0(\x\~ a ) Tailofx P <oo 

Therefore, for any e G (0, 1), there is an I G [L, oo) such that 2 d+a c x p < e for all \x\ > I. 
Then, for \x\ >£, (I2.37P implies 

G p (x) < G p (u) pD{v — u) G p {x — v) 



{u,v} 

(M<§W<M<§|*|) 



+ Y Gp(u)pD{v-u)Gp{ x — V) 



{u,v}dZ d 

{\u\<\\x\, H>|N) 
< c x p sup G p (x — v) + XpP sup D{v — u) 

«:|«|<i|x| {u,v}cZ d 

(\u\<l\x\, \v\>±\x\) 

<2~ d ~ a e sup G p (v) + -^, (2.44) 

u:|u|>||a:| IFIIl 

for some C p = 0(Xp)- If 2£ < \x\ < A£, then we use (12.440 twice to obtain 
G p (x) < {2- d ~ a ef sup G p (v) + 2~ d ~ a e—^—- + 



v:\v\>±\x\ IF 

Cj 

\\d+a 



d+a in llld+Q 

Jb T 



= (2~ d ~ a e) 2 sup G p (v) + (l + e)—^. (2.45) 

v:\v\>±\x\ IF I 

In general, if 2 r I < \x\ < 2 1 for some n G N, then we repeatedly use (12.440 to obtain 
G p (x) < (2~ d - a e) n sup G p (v) + (1 + e ■ • • + e n - 1 " p 



■2" 

Od+a 



|||„|||a!+a 

v:\v\>Jrr\x\ IIFIIlL 



c 

- w^ Xp+ (i- e )iNir- (2 ' 46) 
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For |x| < £, we use the trivial inequality G p (x) < Xp < ^ +a Xp/IMlL +a - This completes 
the proof of f[2"736]) . where K p = £ d+a x P + C p /{1 - e). ■ 



3 Proof of the main result 

In this section, we prove the asymptotic behavior (I1.27P of G Pc in high dimensions. To 
do so, we show in Section [3721 that, if d > d c and L ^> 1, then G p for p < p c obeys the 
same bound as in (12.11) on S p for p < 1. Then, in Section I3~3"l we show that the obtained 
infrared bound on G Pc implies its asymptotic expression (JT727J). The proofs rely on the 
lace expansion (11.121) for G p . 

3.1 Bounds on II p assuming the infrared bound on G p 

In this subsection, we assume the infrared bound on G p and prove bounds on LJ p and 
related quantities, such as its sum IJ p (0) = ^2 x n p (x), in high dimensions. Before stating 
this more precisely, we need introduce the following parameter for a > 0, a ^ 2 and 
d > a A 2 (cf., flUE}): 

A = SU /^ = 0(raA2) ' (3 - 1} 

III x III L 

Proposition 3.1. Let o > 0, a / 2 and d > d c , and assume the properties U.16]) and 
U.18\) in Assumption \l.l[ Suppose that 

P < 3, G p (x) < 3A|||x|||f 2 ~ d [x ^ o}. (3.2) 

If X ^ 1 (i.e., L ^> 1), then, for any x G 

{pD*G p ){x)<0{\)\\x\\r 2 -\ (3.3) 

\n p {x) - 8 0jX \ < 0(L- d )5 , x + 0(\ e )l\xf L aA2 - d)e , (3.4) 

where £ = 2 for percolation and £ = 3 for SAW and the Ising model. As a result, 

n p (0) = l + O(L- d ), (3.5) 



(3.6) 



V aA2 II p (0) = lim 

\k\^o l-D(k) 

[a < 2], 

jzE x \x\ 2 n p (x) = 0(L- d ^) [a>2]. 

We prove this proposition by using the following lemma, which is an improved version 
of [TBI Proposition 1.7]. 

Lemma 3.2. (i) For any a > b > with a + b > d, there is an L-independent constant 
C = C(a, b, d) < oo such that 

Vlllx- y |||-||| y |||- b <|^ aW ^ [a>rf] ' (3 7) 

^ in yiiiL myiiiL - ^11^11^6 [ a< d]. 
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Figure 1: The left figure is an example of the lace-expansion diagrams for percolation, 
and the right one is for SAW and the Ising model. The number I of disjoint paths from 
o to x using different sets of line segments is 2 in the left figure and 3 in the right figure. 



(ii) Let f and g be functions on U 1 , with g being Z d -symmetric. Suppose that there are 
C±, C*2, C3 > and p > such that 

f(x) = d IMIf 2 ~\ \g(x)\ < C 2 5 0>x + C 3 \jx\H d - p . 

Then there is a p' G (0, p A 2) such that, for d > a A 2, 

it \ ( \ \\g i 0(C\C 3 ) . . 

(/ * 9)(x) = ,,, |||d _ aA2 + d - aA 2 +f y ( 3 - 8 ) 

t\ x \h NIL 

Proof of Proposition IJ.il First we note that 

_ om _ om\\x\\i^ o(A) 

1 J ~ ~ lll^lll d - aA2 - III -rill d-«A2- 

III III I, III x III L III x III L 

We also note that the identity G p (y) = 5 0tV + G p (y)t^ j holds for all three models. 
Therefore, by using the assumed bound (13 .2p and Lemma I3.2( i). we obtain (I3.3P as 



(D * G p )(x) = D(x) + Y,D{x- y) G p (y) 



< °^ y 0(L a ) 3A 0(A) 

— IN |||(i-aA2 ' IN _ „,|||d+a |||„,|||d-aA2 — ||| ||| d-c*A2 ' v -- 1 -"; 

ill x ill l y£Z d m x yllli illy III L III X III L 

The inequality H3.4[) is obtained by repeatedly applying (l3.2p -f l3T3|) and Lemma I3.2( i) 
to the diagrammatic bounds on LT p (x) in [18], [26] (II p (x) in this paper equals 5 0tX + U z (x) in 
rrSJ Proposition 1.8]), where I is the number of disjoint paths in the diagrams from o to x 
(cf., Figured]). The proof is quite similar to jTSJ Proposition 1.8] and Proposition 3.1]; 
the only difference is the use of ||| • \\l instead of ||| • |||i and Lemma I3.2( i). Because of this, 
we gain the factor 0(L~ d ) (= 0(A)|||o|||2 A2 ~ d ) in (I3T41) . which is much smaller than 0(A) 



as claimed in [T8l [26] . 

It remains to prove (l3.5l) -( l3T6l) . By (13.41) . we readily obtain (13. 5ft as 

i7 p (0) = n p( x ) = 1 + + O(L-W-V) = 1 + 0(L~ d ). (3.11) 
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Moreover, 

\n p (o) - n p (k)\ 



5^ (i - cosO • x))n p [ 



X) 



< 



o{\") E 



1 — cos(fc • x) 

i(d-QA2)£ 



(3.12) 



If a < 2, then there is a <5 G (0, (2 — a) A ((£ — — d c ))) such that 1 - cos(/c • x) < 
0(\k ■ x\ a+s ), hence 



\n p (o)-n p (k)\<o{\k\ 



IT* £ 

x:|x|<L 
^_^-l) 4 «+*)j A .| a +* < 



a;:|a;|>L 



o+(5 



a;|(«*-o)/ 



(3.13) 



If a > 2, then there is a 5 e (0, 2 A ((£ — — d c ))) such that 1 - cos(/c • x) = \\k ■ x\ 2 + 
0(\k ■ x\ 2+s ) and therefore 



4(0) - n p (k) = \Y,\ k ' x \" n ^) + o{L-")\k\™ E 



\x 



2+8 



x&2 



\x 



\(d-2)e 



2d 



E \x\ 2 n p {x) +0(L' 



-d{l-l)+2+5^]2+5 



(3.14) 



Then, by the above estimates and (11. 18ft . we obtain 

n p (o) - n p (k) = (o(L- d ^+ s )\k\ s [a < 2], 

1 - D(k) ~\jsJ2 x \x\ 2 n p (x) + 0(L-W-V+ s )\k\ s [a>2], 
hence ( 13.61) by taking \k\ — > 0. This completes the proof of Proposition 13.11 



(3.15) 



Proof of Lemma \3.Si The proof of ( 13. 7p is almost identical to that of [TBI Proposition 
l.T(i)] . However, since we are using ||| • |||x, rather than ||| • |||i as in [18], we can gain the 
extra factor L d ~ a for a > d in (13.71) . To clarify this, we include the proof here. First of 
all, since a > b, we have 

E - y\n a Ml" < E iiiz-2/iiriNiz 6 + E i*-vii a Mi b 

y eZ d y-\x-y\<\y\ v-\x-v\>\v\ 

<2 E ™* " vl? Ml- (3-16) 

y.\x-y\<\y\ 

Since \x — y\ < \y\ implies \y\ > we obtain that, for a > d, 

E in* - y\\T Ml ^ 26 Nz 6 El* - y\\T = § Ld_0 iN«z 6 - ( 3 - 17 ) 

y-\x-y\<\y\ y£Z d 



19 



For a < d, on the other hand, we use the identity 1 = l^y^ 3 ^j-j. + l{|,,|>||a;|} an d the fact 
that \y\ > implies \x — y\ > \\y\. Then, we obtain 



\i x - 

v-\x-y\<\y\ 
< 21MIZ 6 



iriMllz 6 



\i x 

V'-\x-y\<\y\ 



\L ai {\y\<l\x\} + 3 a Yl 



—a—b-n 

l l {\y\>l\x\} 



V-\x-y\<\ 



< 2 b \\x\ll b 



E 



\x 



-a + 3 a ^ 



< — \\\x\\ 



(3.18) 



y-\x-y\<%\x\ 



?/:M>§M 



This completes the proof of (13. 7\i . 

The proof of (13.81) is also quite similar to that of [T8J Proposition 1.7(h)], where [18j 
(5.8)] is used. However, [181 (5-8)] is valid only for d > 4, not d > 2 as claimed in [TH 
Proposition 1.7(h)]. In fact, it is not difficult to avoid this problem, and we include the 
proof here to clarify this. First we note that 



(/ * g){x) = \\ghf(x) + 9(v) {f( x -v)- /(*))■ 



(3.19) 



To prove (I3.8p . it suffices to show that the sum in the right-hand side is the error term in 
(13. 8p . For that, we split the sum into the following three sums: 



E = E + E + E -E+E"+E"' < 3 - 20 > 

y:\y\<±\x\ !/:|x-i/|<||x| y:\y\A\x-y\>i\x\ V V V 



It is not difficult to estimate the last two sums, as 

y^'g(y)(f(x-y)-f(x)) < ° {C;i) 
y 



\x 



\d+p 



\\\X\\\ T 



3) 

a/\2+p ' 



y\x-y\<\\x\ 



and 



\d-aA2 



d+p 



d-aA2+p ' 



(3.21) 



(3.22) 



ll|J - |llL y.\y\>l\*\ ll|y|lli 
To estimate the sum ^2 y , we use the Z d -symmetry of g to obtain 

-»)-/(«))= E »(»)( /(J+l0 + /(x ~' -/(*: 



Notice that 

f(x + y) + f(x-y) 



V-0<\y\<i\x\ 



, O(Ci) 11 
< , — ^ x 



|d-oA2 



M 2 /|*| 2 [|x| > |L]. 



(3.23) 
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To verify this for \x\ < §L, we simply bound each / by 0(Ci)|||x 



\aA2- 



For \x\ > \L, 



since |a;±y| > \x\ — \y\ > > L, we have f(x±y) = C\\x ±y\ . Then, by Taylor's 
theorem, since I ± 2^ + I < ? < 1, we have 



\x±y\ aA2 - d = \x\ aA2 ~ d \ 1±2 



x-y , \y\ 



4\ (aA2-<2)/2 



Ixl 2 Ixl 4 



Lr 



r I 12 

-,\ x ' y i q' \y 



1 =F (d — a A 2)- 



.r 



(3.24) 



and (I3.23P follows. Therefore, if |x| < |L, then |?/| < ~L and we obtain 



£'<Kl/)(/(* -*)-/(*)) 



< 



0{C 



\\d-aA2 



E 



III- — * < L d+ P 



C 3 < 0(CiC 3 ) 



|d-aA2+p ' 



(3.25) 



If |a;| > |L, then \x\l = \x\ and we obtain 



^g(y){f(x-y)-f(x)) 



< 



0{C 



\d-aA2+2 



\\-f 



E iv 



„.i2/ Caljjj^} C 3 l{|y|>L} 



Id+p 



j/:0<|i/|<i|a!| 

L-" +2 [p>2], 

:^Slxilog|x| [p = 2], 

|a;| 2 ^ [p<2]. 



x 



|d-aA2+2 



(3.26) 



Summarizing the above yields (13. 8p . This completes the proof of Lemma [3.21 



3.2 Proof of the infrared bound on G p 

In this subsection, we prove that the hypothesis of Proposition 13. II indeed holds for p < p c 
in high dimensions. The precise statement is the following: 



Theorem 3.3. Let a > 0, a 2 and d > d c , and assume the properties U.16\) . U.18\) 
and PJ%) in Assumption \l.li Then, for L ^> 1 and p < p c , 



G p {x) < 0(L- 



aA 2) I,, ,|| *A2-d 



[x ^ o] 



(3.27) 



Proof. Let 



g p = pV sup 



G p (x) 



/ \ III -T- III aA2-d ' 
x t z ° '* HI Jb II r 

where we recall the definition (13. ip of A. Suppose that the following properties hold: 

(i) g p is continuous (and nondecreasing) in p £ [l,p c )- 

(ii) <?i < 1. 



(3.28) 
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(iii) If A 1 (i.e., L>1), then g p < 3 implies g p < 2 for every p G (l,p c )- 

If the above properties hold, then in fact g p < 2 for all p < p c , as long as d > d c and 
A < 1. In particular, G p (x) < 2A|||x|||2 A2_d for all x ^ o and p <p c (< 2). By Lemma [221 
we can extend this bound up to p = p c , hence the proof completed. 
Now we verify those properties (i)-(iii). 

Verification of (i). It suffices to show that, for every p G (l,p c ), sup^ G p {x) /\x\°^ 2 ~ d 
is continuous in p G [l,Po]- By the monotonicity of G p (x) in p < po and using Lemma [2~4| 
we have 



G p {x) ^ G P0 {x) ^ K P0 \\x\\- L d - a _ K 



-d—a 

Ml _||| OS/V2 — ci — ||i III q:A2— rf — hi ||iaA2-d ||| _ ||| a+aA2 ' \<J-^o) 

III III L III x III L III x III L III x III L 

On the other hand, for any xo ^ o with D(xq) > 0, there exists an R = R(p ,xo) < oo 
such that, for all \x\ > R, 

|a+aA2 - iii laA2-d- (3.30) 

\l IFoIl 



\\x 



Moreover, by using p > 1 and the lower bound of the second inequality in (12.251) . we have 

QA2-d — iii III Q!A2— d — ni |||aA2-d" ^O.Olj 



iFoli IIFoIl I x o||Il 

As a result, for any p G [l,£>o]> we obtain 

G p (x) G p (xo) G p (x) 
sup AO , = AO , V max .„ , . 3.32) 

x /o |||x|||£ A fflxolir v-°<\*\<R \lx\l a L A2 - d V ; 

Since G p (x) is continuous in p (cf., Lemma f2.2[) and the maximum of finitely many con- 
tinuous functions is continuous, we can conclude that g p is continuous in p G [l,Po]> as 
required. 

Verification of (ii). By the first inequality in (I2.25j) and the definition (13. ip of A, we 
readily obtain 



9i = 1 V sup , IH _2 A i_ d < 1 V sup ')J 2 _ d = 1, (3.33) 



GAx) SAx) 

\t . < 1 V sup ^- ' 

x^o A HI X |||^ x^to /V HI X |||£ 

as required. 



Verification of (iii). li d > d c , X <^ 1 and g p < 3, then, by Proposition 13.11 H p satisfies 
(13.4p - fl3.6p as well as (13.151) . We use these estimates and the lace expansion to prove 
g p < 2 as follows. 

First we recall ( 02]) and (OTj) : 

G p = 17 p + II p *pD* G p , S p = 5 + pD* S p , (3.34) 
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or equivalently 

n p = G p *(5-n p * P D), 6 =(5- pD) * S p . (3.35) 

Inspired by the similarity of the above identities, we approximate G p to rJJ p * S q with 
some constant r G (0, oo) and the parameter change q G [0, 1]. Rewrite G p as follows: 

G p = rll p * S q + G p * 5 — rlJp * S q 

= rllp * S q + G p * (5 - qD) * S q - rG p * (5 - n p * pD) * S q 

= rI7 p *S q + G p * E p ^ r * S q , (3.36) 

where 

Ep,^ = 5 - qD -r(5 - n p * pD). (3.37) 

We choose q, r to satisfy 

E P>q A°) = V QA2 £ W (0) = 0, (3.38) 

or equivalently 

l-q-r(l-n p (0)p) = 0, 
-q + r(n p (0) + V aA2 n p (0))p = 0. 

Solving these simultaneous equations for r and using (13.61) . we obtain 



(3.39) 



r = (i+ P v aA2 n p (o)y l = 1 + 1^ 



[a < 2], 

0(L- d ^) [a>2]. 



(3.40) 



On the other hand, by taking the Fourier transform of (13.341) and setting k = 0, we obtain 

x P = n P (o) + n p (o)p Xp , (3.4i) 

or equivalently IJ p (0) = +PXp)> an d therefore 

q = i _ r (i _ n p (0)p) = 1 - G (0, 1], (3.42) 

J- + PXp 

where we have used p > 1, x P > 1 and (13.401) to guarantee the positivity (by taking L 3> 1 
if a > 2). 

In addition, by solving (I3.4ip for \v an d using (13.51) . we have 

4(Q) = i + o(^) 

1-/7,(0)? l-/7 p (0)p 

hence 1 — n p (0)p > 0. In particular, p < iJ p (0) _1 = 1 + 0(L~ d ) < 2, as required. 

It remains to prove G p (x) < 2A|||x|||2 A2_d . To do so, we use the following property of 

Ep,q,r'- 
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Proposition 3.4. Let q,r be defined as in ( KJW - [^J^ ). Under the hypothesis of Propo- 
sition IKT\ there is a p G (0, a A 2) such that 



\{E m>r * S q ){x)\ < 0(L~ d ^)(l {a>2} S , x + 



\d+P / " 



(3.44) 



For now we assume this proposition and complete verifying the property (iii). First, 
by rearranging ( I3.36|) and using S q < S\ as well as ( 13. 5 P and (I3.40p for L ^> 1, we obtain 



G p = rLJp * S q + G p * Ep^ r * S q 

= rLf p (0)S q - r(n p (0)5 - LI p ) *S q + G p * E Ptq>r * S q 

< (1 + 0(L- d ))S 1 - r(n p (0)5 - R p ) *S q + G p * E p>q , * S q . 

Then, by Proposition 13.41 and Lemma I3.2( i). the third term is bounded as 

IP 



\(G P * E p , q , * S q )(x)\ < 0(L~ d ^) £ + 



i \w d +p 
\ x -y\h 



0(L' d ^)X 



< 



\\x 



\d-aA2 



Also, by (13 .4p and Lemma l3~27 i). the second term in (I3.45P is bounded as 



((n p (o)6 - Lip) * s q ){x)\ = ^2nM(s q (x) - s q (x - y )) 

< \n P (y)\ s q (x) + \ n p(y)\ s «( x - y) 

y+o y+o 
0(L^~ 1 ))A 



< 



\x 



\d-a/\2 



Putting these estimates back into (13.451) . we obtain that, for L> 1 

A 



G p (x) < (l + 0(L- d ))- 



\x 



\d-aA2 



0(L~ d ^)\ < 2A 



x 



\d-aA2 



\X\ 



d-a/\2 ' 



(3.45) 



(3.46) 



(3.47) 



(3.48) 



as required. This completes the proof of Theorem 13.31 assuming Proposition 13.41 



Proof of Proposition \3J\ First, by substituting q — 1 — r(l — LI p (0)p) (cf., ( 13.421) ) into 
(l3T37j) and using 1 - r = prV aA2 n p {0) (cf., (jg^Djl ). we obtain 

Ep, q , r = pr(v aA2 n p (0)(5 - D) - (n p (0)5 - Lip) * £>). (3.49) 

Using this representation, we prove (I3.44p for \x\ < 2L and \x\ > 2L, separately. 
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For \x\ < 2L, we simply use ( 12.1 p to bound \(E Ptq<r * S q )(x)\ by 

\{E M>r * S q )(x)\ < \E p>q>r (x)\ + 0(L~ d ) \E P , q ,r{y)\- 



(3.50) 



By (I3T49D . we have 

I Ep q r (x 



pr 



< \v aA2 n p (0)\(s o , x + d(x)) + (n p (o)6 -n p )*D 



x 



= I V QA2 i7 p (0) | (S , x + D(x)) + Yl n p( z ) i D ( x ) ~ D ( x ~ z )) 

z^o 

< I V QA2 i7 p (0) | (6 0>x + D(x)) + \ n p( z ) I ( D ( x ) + ^ - *)) • 

Using (|3Ti|l - (j3T6p and igM) , we obtain that 

\E p>q , r (x)\ < 0(L- d ^)t {a>2} (5 0tX + 0(L- d ))+ 0(L- d )J2\ n p( z )\ 

z^o 

<0(L- d ^)l {a>2} S , x + 0(L- M ), 
and that, by summing (13.5 1[) over x G Z d , 

£ l*W(z)| < o(L- d )(2\v^ 2 n p (o)\ +2j2\n P (z)\) < o(l- 

x& d ^ z^o ' 

Therefore, for |x| < 2L, 

\(E m , r * S q )(x)\ < 0(L-W-V)l {a>2} 5 0>x + 0(L~ 



(3.51) 



dl X 



< 0(L-W-V)l {a>2} 5 0>x 



0{L 



-d(e-i)+ P \ 



\x 



\d+P 



(3.52) 



It remains to prove (I3.44p for \x\> 2L. To do so, we first rewrite (E p ^ r * S q )(x) as 

ld]„ „—ik-x 



(Ep, q ,r * S q )(x) 



Ep, q , r {k) 



(3.53) 



Then we split the integral with respect to t into J Q T and J^°, where T is arbitrary for now, 
but it will be determined shortly. For the latter integral, we use the Fourier transform of 
( E3SD , which is 



E m , r (k) =pr(l- D(k)) ( V- 2 77 p (0) - ^l^ m ) . 



(3.54) 
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Because of (11.181) . (13. 6p and (I3.15p . there is a 5 > such that 

E Mtr (k) = 0(L~ d{£ ^ +aA2+5 )\k\ aA2+s [a ^ 2]. 



(3.55) 



Since 1 — qD[k) > q(l — D(k)), the contribution to (13.531) from the large-t integral is 
bounded as 



dt 



< 0(L- 



d d k 



E P,gA k ) e 



-t(l-qD{k))-ik-x 



-l)+aA2+<5\ 



dt 



d d k 

[-7r,7T] d ( 27T )° 



1*1 



aA2+S -tq(l-D(k)) 



(3.56) 



Since p > 1, we have g > 1 — r/(l + xi) > 1 — r/2 (cf., (I3.42p ). which is bounded away 
from zero when L ^> 1. Therefore, by using (I1.18p . we obtain 



d d k 



aA2.\-l- 



\k\<*™+> e -mi-»w) = o(L t) 



d+S 



(3.57) 



hence 



dt 



d d k 



T,7T]d ( 2 ^ 



-t(l-qt)(k))-ik-x 



<0{L- M )T-^, 



(3.58) 



Let 



[a A 2)8 



d + af\2 + 5 



aA2~p 



(3.59) 



Then, since Ixl > 2L, 



\x 



\d+P 



(3.60) 



To estimate the contribution to (I3.53P from the small-t integral, we use the identity 



dt 



d d k 



-7r,7r] d (2tT 



E p>gr (k) e ~t(i-qD(k))-ik-x 



n=0 



^-^ nl 

(3.61) 



where, by ( ESD and f l3~6|) . 



(S M , r *D* n )(z)= prV QA2 77 p (0) 

" v ' d 

0(L-^-D)l {Q>2} 

- pr R p(y) {D* {n+1) (x) - D*( n+ ^ (x-y)). (3.62) 
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In the following, we use the decomposition (I3.20p of J2 y an d estimate the contribution to 
( 13.611) from Yl' y , J2'y and Ylyi separately. 

First we estimate the contribution from J2 y = lt2y\x-y\<±\x\- Since \y\ > \x\ — \x — y\ > 



||a;| in this domain of summation, we bound \IJ p (y)\ by 0(A^)|||x|||^ QA2 (cf., ( 13.41) ) and 
then use (Olj) . £j 1 < OQ|:r|£) and Y! y D< n+l \x - y) < 1. As a result, 

E*^)(^'>(») - - »)) < + 1) 

oi X 7- \ III III L / 



< 



0(L- d ^ +aA2 ) (0{L 



aA2^ 



\\.r 



\\d+af\2 



Lr 



||aA2 



Similarly, for a > 2, 



0(L 



< 



< 



0(L-^- 1 )+ Q ) fO[L 2 )n 



\\d+a 



V III x II 



||d+2 



v 



/i 



.r 



(3.63) 



+ 11. (3.64) 



To estimate the contribution to (13.611) from J2y = Z)j, : |y|A|x-i/|>i|x| 



in (13.62}) , we bound 

D< n+1 \x) and D'^fc-j/) by 0(L aA2 )n/\\x\\ d L +aA2 and then use (S3} to bound |/7 p (i/)|. 
The result is 

Yl' ' n p {y){D< n+1 \x) - D< n+1 \x - y)) 



< 



0(L aA2 )n 



Lr 



||d+»A2 



£ 1^(^)1 



< 



Lr 



L v:|»l>i|x| 

0(L -(£-l)(aA2) )rl 
|id+2(aA2)+(f-l)(d-d c ) 



< 



£U£-d(^-l)+2(aA2)) 



n 



Lr 



iid+2(aA2) 



(3.65) 



Similarly, for a > 2, 



0(1/ 



-d(£-l)- 



Y^"D(y)(D-(x)-D'"{x-y)) 



Lr 



id+2 



< 



0(L-^" 1 )+ 4 ) 



n 



Li" 



(3.66) 



Finally we estimate the contribution to (13.61 j) from = Yl y -\y\<±\ x \ ^ n (EB2J). By the 
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Z d -symmetry of II p and using (13. 4p and the assumption (11.241) . we obtain 
Y^n p {y){D* [n+l \x) - D< n+l \x - y)) 



\d+aA2+2 / j 



11' 



l/:|j/l<*W 



(d-aA2)£ 



< 



< 



0(^-(/-l)(aA2)) n 



id+2-(d-aA2)£ 



\\x 



\d+af\2+2 



[d + 2 > (d - a A 2)1] 
[d + 2 = (d - a A 2)£] 
[d + 2 < (d — a A 2)1] 



T j d+2-(d-aA2)l 



0(L -d(£-l) + 2(aA2) )n 



.r 



ia!+2(aA2) 



(3.67) 



where, to obtain the last inequality for d + 2= (d — a A 2)£, which implies a < 2, we have 
used fact that (|||x|||i/L) Q: ~ 2 (l + log(|||x|||i/L)) is bounded. Similarly, for a > 2, 



O(L-^-D) 



0(L" 



d(^-l) 



D(y)(D"'(x) - D"'( X - y)) 



< 



0(L- d ^ +2 )n 



\\x 



|d+4 



l/:|l/|<*M 



0(U 



0(1 



-d(e-l)+4 



) n 



d+a 
L 



|(i+4 



(3.68) 



0(L2) 



Now, by putting these estimates back into (13.621) . we obtain 



\(E M , r *D*»)(x)\< 



( L -d(/-l)+oA2) / 0(L 



qA2> 



I <2+aA2 



|aA2 



-n + 1 



(3.69) 



hence, by (I3.59p . 



n=0 



< 



0(1 



-d(l-l)+aA2t 



X 



|d+aA2 



0(1 



aA2\ 



\X 



\aA2 



-T 2 + T 



0(I~ d ^^ +ah2 ) ^ 



|d+«A2 



0(1 



-d(t-l)+p\ 



\x 



\d+p 



(3.70) 



This completes the proof of Proposition [3T 
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3.3 Derivation of the asymptotics of G Pc 

Finally we derive the asymptotic expression fjl.27J) for G Pc . First, by repeatedly applying 
(13. 36ft . we obtain 



G p — rllp * S q + G p * E p q r * S q 

rll p * S q + (rlJp * Sq + Gp* E p ^ r * S q ) * E p ^ r * S q 
rn v *S q *(5 + Ep^ r * Sq) + G p * (E M>r * S q )* 2 



N-l 



rllp * Sq * * S q )* n + G p * (E p ^ r * S q )* N . (3.71) 



n=0 



By Proposition 13.41 and Lemma l3.2( i). we have that, for p < p c , 

\(E p , q , r * S q y n (x)\ < 0(L~ d ^ n ) (l {a>2} 8 0>x + ) , 



hence, for any iVeN, 



JV-l 



E l^w * S q )™{x)\ < (1 + 0{L- d ^))5 , x + ° {L , +p " ; . (3.73) 

n=0 IIfIIIl 

Therefore, we can take iV — > oo to obtain that, for p < p c , 

oo 

Gp = rllp * Sq * ^(E M , r * S q )* n = Hp * Sq, (3.74) 

n=0 

where, by ([33} and fl3T73D . 

H p (x) =r(n p * f2{E p , q , r * S g )*A (x) 

^ n=0 ' 



-d(€-l)+p>l 



x ( (1 + 0(L"^-)))^ + ) • (3.75) 



- y\\ d L +p 

Notice that, by Lemma l3.2( i) and using (13.401) and d + p<(d — aA 2)£, 

H p (x) = (r + 0(L- d ))5 , x + U[ d+p } . (3.76) 

NIL 

Now we set p = p c , so, by (13. 42ft . q = 1. By Proposition 12.11 and Lemma l3T27 ii). we 
obtain the asymptotic expression 

IIP Hi l x l \\\ x \\\l 
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Since H Pc is absolutely summable, we can change the order of the limit and the sum as 

oo 

H Pc (0) = lim H Pc (k) = lim rf[ pc (k)J2 {E Pc ,iAk) Sx(k)) n 

71=0 



OO 

n Pc (o) + rn Pc (o)J2( ,ir 4ci,r(fc)^i(A;)) n (3.78) 



n=l 



By (I3.43P and the fact that Xp diverges as p •f p c , we have IJ Pc (0) — p c 1 . Moreover, by 
fl334j) and fl3~T5|) . 



4 cjl , r (A;)^(A;)=4 ,i,r(A;)(l-D(A;)) 1 ->• 0. (3.79) 
Therefore, 

A = H^O)- 1 = ^ = p c (l +p c V aA2 i7 Pc (0)). (3.80) 

This completes the proof of Theorem 11.21 ■ 



A Verification of Assumption 11.1 



In this appendix, we show that the Z d -symmetric 1-step distribution D in f ll.25p . defined 
more precisely below, satisfies the properties ( 11. 16j) . ( 11.181) . f ll.20p . ( 11.21 jl and (1 1.24ft in 
Assumption 11.11 

First, for a > and a ^ 2, we define 

T a (t)= * Xa 'l [teN], (A.l) 

Z^seN * 

Next, let h be a nonnegative bounded function on M. d that is piecewise continuous, Z d - 
symmetric, supported in [—1, l] d and normalized (i.e., J^ 11 ^ d h(x)d d x = 1); e.g., = 

2 _<i l{|| ;I ,j| oo < 1 }.. Then, for large L (to ensure positivity of the denominator), we define 

Ul{x) = hix/ h ^ /n [x e Z% (A.2) 

E ye z<* Ky/ L ) 

where (cf., [121 122]) 

<r£ = £ |x| 2 f/ L (x) = 0(L 2 ), (A.3) 

€(-l + A,l-A) m>'l\ 

= i-Si*|» + o((i|*|)«<) [| fcH0 |, (A - 4) 
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for some A £ (0, 1) and £ £ (0, 2). Combining these distributions, we define D as 



teN 



We note that the above definition is a discrete version of the transition kernel for 
the so-called subordinate process (e.g., [8]). Just like ( 1A.5[) . the transition kernel for 



the subordinate process is given by an integral of the Gaussian density with respect to 
the 1-dimensional a/2-stable distribution. Bogdan and Jakubowski j9] make the most 
of this integral representation to estimate derivatives of the transition kernel. This is 
close to what we want: to prove (11.2411 . However, under the current discrete space-time 
setting, we cannot simply adopt their proof to show (11.24)1 . To overcome this difficulty, we 
will approximate the lattice distribution in (1A.5|1 to a Gaussian density (multiplied 



by a polynomial) by using a discrete version of the Cramer-Edgeworth expansion [TJ 
Corollary 22.3]. 

Before doing so, we first show that the above D satisfies ( 11. 18ft and f l 1 . 201) . 

Verification of $1.180 and ( f 1.200 . Due to the above definition of Ul, we can follow the 
same argument as in [221 Appendix A] to verify the bound on 1 — D in ( 11.201) . Moreover, 
if (11.1811 is also verified, then we can follow the same argument as in [TUl Appendix A] to 
confirm the bound on HD*"^ in (1 1 . 2 [) as well. 

It remains to verify (I1.18P for small k. First we note that 

t 

1 - D(k) = £(1 - U*) T a (t) = (1 - U) T a (t) £ U s ~\ (A.6) 
teN teN s=i 

where U is an abbreviation for U^k). If a > 2, we can take any £ £ (0, a/2 — 1) to obtain 



i - D(k) = (i - u) T a (t) 1 - (i - u) E D 1 - 

teN s=i teN s=i 

= (l-U) y £tT a (t) + 0{(l-U) 1+ t) 1 (A.7) 



teN 

where we have used the inequality 

teN s=i teN s=i ^ U ' 

< 2^(1 - ?7)^t 1+ «T Q (t) = 0((1 - f>) 5 ). (A.8) 

teN 

This together with ( IA"l3]l -( EOj) implies < TT7TH|> for a > 2, with e = ( A (2$) and 



teN 
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If a G (0, 2), on the other hand, we first rewrite (1A.6j) for small k by setting u = log 1/U 
and changing the order of summations as 

i-dw^I/E^E^^ V v ^ — • (A - 10) 

We note that, for small k, 

LzE = l-U + 0((l-U) 2 ), u = 1-U + 0((1-U) 2 ). (A.ll) 

Therefore, by a Riemann-sum approximation, we can estimate the numerator in (1A.10I) 

as 

oo / ^ -l-a/2 



(*>«) 



/OO /*OC 
ds e" s y dt t" 1 ""/ 2 

-r(l - a/2) u a/2 ^ (1 + 0(u)). (A.12) 



This together with flA~3l)-flA~4l) and flA~9l)- flA~T2ll implies f lTTSD for a G (0, 2), with e = £ 
and 

2 ^rW). (A.13) 



This verifies that D in f TA~5|) satisfies both ffTTT8|) and f OUj) . 



Verification of U.16\) . U.21\) and \1.2J$ . To verify these x-space bounds on the transi- 
tion probability D* n and its discrete derivative, we use the Cramer-Edgeworth expansion 
to approximate the lattice distribution U^(x) in flA.5[) to the Gaussian density u a 2 t (x) 

(multiplied by a polynomial of xj \J a\i) , where 



u c (x)= [J-Y^f-^L). (A.14) 



d/2 

2-kc J ^ V 2c 

Before showing a precise statement (cf., Theorem IA.1I below), we explain heuristic 
derivation of the formal expansion (1A.21j) . First, we recall the expansion in terms of 
cumulants Qa for n G Z+: 

_ jL (jU \n s 

\o g u L {k)= ^W^r- (A.i5) 

(ll»lll>l) 
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Since Ul is Z d -symmetric, we have = if ||n||i is odd, and Q(2,o,...,o) = ■ ■ ■ = Q(o,...,o,2) = 
(J 2 L /d. Therefore, 

2 00 d / ■ i \ n 

iog^) = -^r+£ e ^n^- (a.iq) 

1=4 Hezf s=1 s ' 

(\\n\\l=l) 

By the Fourier inversion, we may rewrite U^(x) as 
U ?(x) = [ 7^ U L (kfe- ik - x 



(ll*s||i=0 

(alt)-^ ( ^- e--l fc l 2 -^exp ( fY^QKi*)), (A.17) 



where, in the third equality, we have replaced fc by k/ \fo\i and used the abbreviations 

<},(«)= E $n^. (A.18) 

(ll»lli=0 



.r 



^ n6 



Notice that the coefficients Qh/^l f° r ll^lli = ^ are uniformly bounded in L. Then, the 
exponential factor involving higher-order cumulants in flA. 17|) may be expanded as 

(00 \ 00 m 

E*" ,/a ^( i *)) = 1 +E^ E n( rir/2 ^ +2 w) 
Z=2 ' m=l ' I lr ..,l m >2r=l 

00 b'/ 2 J m 

= l + E rJ/2 E^T E U^ + 2(ik). (A.19) 



Let 



771! 

j=2 m=l Ji,— ,«m>2 r=l 

(/l + -+im=j) 



Li/2J 



PoW = l, P 1 (ifc) = 0, P^A;) = ^ ^ ^ f[Qi r+ 2{ik) [j > 2]. 



771! 

771=1 il,---j'm>2 r=l 

(«l+-+«m=j) 



(A.20) 

Then, by (1A.17P and (1A. 19j) . we arrive at the formal Cramer-Edgeworth expansion 

_ 7^ e-^ |fc|a -*' fi E ri/2p i( <A! )- (A. 21) 
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Now we note that, if \fo^t\— ir, ir} d is replaced by M d , if ^°1 * s re pl ace d by X^=o ^ or 
some i < oo, and if x is considered to be an element of ~R d instead of Z d , then we obtain 



d d k 



d d k 



(2ir) d 
y d/2 y j/2 p f 



(2tt) £ 



e 2d\ k \ 2 



ik-x 



(A.22) 



3=0 



where P 3 is the differential operator defined by replacing each ik s of Pj(ik) in (1A.20[) by 

— d/dx s : 



P = 1, Pi = 0, 

Notice that, by f lAl8|) and il) . 



-9 



-9 



P. - p 



[J > 2]. 



(A.23) 



[olty^PjV^i) in ( 1A~22]) = # 



r2i 



(A.24) 



where -ff|+ 2 i s a polynomial of degree at least j + 2 and at most 2j (due to the symmetry 
of U£). The coefficients of the polynomial are uniformly bounded in L, as explained below 
flA~T8D . 

The following theorem is a version of Corollary 22.3] for symmetric distributions, 
which gives a bound on the difference between U^(x) and (lA.22j) . 



Theorem A.l. For any x E 7h d , t G N and £ G Z 



+ • 



(1 + \x 



j=0 



< 0(L" d )t 



c(N|,-(d+£)/2 



(A.25) 



where x and Pj are defined in liA.18\) and jiA.2S\) . respectively. 



Before using this theorem to verify (11.161) , (11.211) and (11.241) , we briefly explain how to 
prove that the contribution which comes from 1 on the left-hand side of (1A.25|) is bounded 
by 0(L~ d )t~^ d+e ^ 2 , as in First, we split the domain of integration in the Fourier 



4 To investigate the effect of the factor |ai| £+2 on the left-hand side of (|A.25j) . we also use identities 
such as 



d d k 



Ql + 2 



d(ik 1 y+ 2 



Ul 



k 



toft 



which is a result of integration by parts. 
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space into E 1 = {k e R d : |fc| < v 7 ^}, ^2 = H d \ ^1 and ^3 = Rd \ E i- Then > 

the difference between U^{x) and (1A.22|) is equal to T\ + X 2 — X 3 , where 



X, = (aft)^ 2 
(^lt)" d/2 



X 



d d fc 

d d k 



-ik-x 



Ui 



3=0 / 



-ik-x 



u, 



x 



(A.26) 
(A.27) 
(A.28) 



Since ( 1A.25I) for t — 1 is trivial, we can assume t > 2 with no loss of generality. Then, it 
is not difficult to prove that X 2 and X3 are both bounded by 0(L~ d )t~ ( - d+e ' > / 2 , due to direct 
computation for X 3 , and due to ( 1A.4j) and similar computation to [TOj (A. 2)] for X 2 . For 



Xi, we can bound the integrand by Ct~ e / 2 (\k\ e+2 + |£;| 2£ )e~ c ' fc ' 2 for some L-independent 
constants C, c G (0, 00), due to a version of [TJ Theorem 9.12] for symmetric distributions. 
Then, by direct computation, we can prove that Xi is also bounded by 0(L~ d )t~( d+e ^ 2 . 

Now we apply (IA.25I) to verify the x-space bounds (jl.lfip . ( I1.2ip and ( ll.24p . In par- 
ticular, by fTA3|) and flAT23|) ffAT2%]) . 

00 £ 

3- 



oo oo £ / 

t=l t=l 7=2 ^ 



E 

i=l 



0{L~ d ) 

t (d+l)/2 



1 A 



\x\ 



T a (t). 



(A.29) 



The leading term is bounded as 

00 

Y,^lt{x)T a {t)<0{L- a 



eX P ( - Si) 



£l+(d+a)/2 



< 0(L" 



^ f-l-(*K»)/2 
l<*<|as/*Till? *>l*/ff£|? 



E 

l<t<l*MJ? 



^l+(d+a)/2 



+ 0(|||X/ ( T L |||7^ ) ) 



0(L Q )|||x|| 



-d—a 



iil • (A.30) 

The second term on the right-hand side of (1A.29|) is bounded, due to ( 1A.24I) . as follows: 
for any j 6 {2, . . . ,£} and h 6 {j + 2, . . . , 2j}, 

h 



E'- 



i/2 



< 0(L 



< 0(L 



-d-h\ 



-d-h\ 



X 



exp f _ 

E V 2 4* ; 1 /-l-(d+/i+j+a)/2 
t l+(d+h+j+a)/2 "'" 
l<t<|||x/a £ |||f t>|*/«ri,l? 

/l O(L d+ ^' +a ) 0(L J+Q ) 0(L 2+Q ) 



id+/i+j+a 



id+j+o 



< 



.r 



|d+2+a ' 



(A.31) 
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Therefore, 

oo I 



EE^'/W * ) -V»T.M < £S (A.32) 
*-i Vv¥/ III x III r, 



4=1 j= 2 W U L 1 / III'*' III i 

Similarly, the third term on the right-hand side of (1A.29|) is bounded as 



t (d+£)/2 I I U 

i=l \ v 1 

<0(L- d ) ((fr)" +2 ^ r(^)/ 2 + ^ r l-(d+<+a)/2 

V^ 1 ^ l<t<|*/<Ti|? t>|«/<T £ |f 

'0(L- d+£+2 )|||x|||^- 2 [d + a>2], 

0(Z,-*+^+ 2 ) HarJI^- 2 log lar/o-ilh [d + a = 2], (A.33) 
O(L a+ 0|||x|||^- a - £ [d+a<2], 

which is further bounded by 0(L a+2 )|||x|||^ d ~ Q:_2 for sufficiently large I. Summarizing the 
above estimates, we can conclude (11.161) : 

D(x) = £ Vait {x) T a (t) + < -±>. (A.34) 

t=i IF Hi III x III l 

The bound ( ll.2ip on the n-step transition probability is then automatically verified, 
due to the argument below (11.211) . Heuristically, since 

oo 

D* n (x) = J2ui t (x)T: n (t), (A.35) 

t=n 

this suggests that 

T* n (t) < 0(n) T aA2 (t). (A.36) 

In fact, we can verify this (or a stronger version) by following the same argument as given 
below ( ll.2ip . but we omit the details here. 

Finally we verify (jl.24p by using ( 1A.25|) with sufficiently large t and ( ]A.34[) — (1A.36|) . 
For \y\ < (so that \x ± y\ > ||x|), we obtain 

™ v _ D m (x + y) + D m (x-y) 
{X> 2 

^ ( 7l\ d/2 ( ,,9, e -v\x+y\ 2 /t : -r,\x-y\ 2 /t\ Q( L aA2+2\ 

where we have set r] = d/{2<j\) = 0(L~ 2 ) for convenience. By a Taylor expansion, 

p-r?|x+j/| 2 /< _L p -»7|a;-s/| 2 A 0(n\v\ 2 ) , ,9, 

2 t ' v ' ; 
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Using this and (1A.36|) and following the same analysis as in flA.30]) — (IA.31j) , we can bound 
the sum in (1A.37I) by 

~ e -n\A 2 /t 0(n- (aA2 ^ 2 )\v\ 2 
Ofa W/2 )li/l 2 "Ei4T^= n - < A ' 39 > 

t=l lll X IHl/v^ 

This together with f lA~37j) and \y\ L = \y\ V L yields ( Oij) . ■ 
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